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Abstract 

For the q— state Potts model new order parameters projecting on a 
group of spins instead of a single spin are introduced. On a Cayley tree 
this allows the physical interpretation of the Potts model at noninteger 
values go of the number of states. The model can be solved recursively. 
This recursion exhibits chaotic behaviour changing qualitatively at critical 
values of go • Using an additional order parameter belonging to a group of 
zero extrapolated size the additional ordering is related to a percolation 
problem. This percolation distinguishes different phases and explains the 
critical indices of percolation class occuring at the Peierls temperature. 
Keywords: Potts model; Cayley tree; Bethe lattice; Logistic equation; 
Percolation; Deterministic Chaos. 

1 Introduction 

An interesting feature of the q-state Potts model ^] is the possibility to con- 
sider its thermodynamic observables as analytic functions of the number q of 
spin states. Extrapolating q to unphysical values (<7 = 0(||)or(7=l/2(j4|]))or 
taking the derivative at q — 1 (||) the observables of other models (resistor net- 
work resp. dilute spin glass or percolation) are obtained. Since this equivalence 
holds only for the thermodynamic observables the connection to the dynamics 
of the original Potts model is lost. Extrapolation requires explicit expressions 
of the observables as functions of q. This can be achieved, if we consider the 
q— state Potts model on a Cayley tree, which can be solved with a recursion 
formalism |6|-]T^|. Extrapolation in q can be done in the recursion for the order 
parameter to obtain relations to new models. As an example we mention the 
result H] that the percolation case q = 1 is equivalent to the logistic equation, 
which describes the formation of percolation clusters. To interpret non-integer 
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values of q we introduce a new order parameter. Instead of distinguishing one of 
the q values of the spin we assume ordering for a whole group of qi spins. A next 
neighbour interaction favouring equal spins will also favour spins in the same 
group. Therefore it is not surprising that we find the same mean field equations 
as in the case q% = 1, however with an effective number of spins q = q/qi and 
a redefined coupling K (K ) as function of the original next neighbour coupling 
Kq. Since q and q\ may be chosen large this allows the interpretation of the 
q— state Potts model with noninteger q in terms of a physical model. Due to 
the redefinition of the coupling not all values of K may correspond to physical 
values of the original coupling Kq. However, an extrapolation in the coupling is 
much easier to perform than the one in the number of components. An example 
will be given. At qo — 1 we encounter the logistic equation exhibiting chaotic 
behaviour depending on the value of K as shown in Q . One of the aims of the 
present paper is to study the transition into chaos for real values of qo. 
A problem on a Cayley tree occurs if the recursion allows for more than one 
stable solution. On a normal lattice surface effects can be neglected and the 
homogeneous system can be characterized by a free energy per site /(to) as 
function of the magnetization. The fixed point solutions of the recursion corre- 
spond to different minima of /. The phase with the lowest value of / will be 
stable, the others are thermodynamically metastable. On a Cayley tree / can- 
not be obtained without further assumptions due to the presence of the border 
points. At a fixed point the interior of a large Cayley tree exhibits the same 
magnetization for all sites. As first treated in Ji], for q = 2 and later on 
generalized for q > 1 by Peruggi |j| a free energy per site / can be found by 
integration. This / describes the interior of a large Cayley tree and can be 
taken as representative for the infinite Bethe lattice. At zero magnetic field one 
obtains II from the absolute minimum of / for ferromagnetic coupling one tran- 
sition from the disordered state into a magnetized state and one transition for 
negative coupling into an antiferromagnetic state. The transition obtained by 
the presumably exact Bethe-Peierls approximation corresponds to a spin- 
odal point of a metastable phase. The picture of a homogeneous free energy 
seems not to be very natural on a Cayley tree since most of the sites are ei- 
ther boundary points or lie in the transition region. As shown in |lJJ and ju| 
the values of the surface magnetizations decide which of the stable fixed point 
is adopted. In this picture the previously as metastable discarded phase ap- 
pearing at the Bethe-Peierls temperature corresponds to a physical phase. To 
reconcile the recursion on a Cayley tree with the approach of a homogeneous 
/ various proposals have been made. The trivial solution [Til by chosing the 
boundary conditions such that / holds also for the Cayley tree does not look 
very natural. In / ought to be computed by the differenc of free energies 
between the tree and a slightly modified lattice which is justified only if the 
solution is homogeneous. Monroe |TJ]] conjectured that the most stable fixed 
point corresponds to the physical phase. This criterion reproduces at zero field 
the phase diagram of Peruggi M . The merit of this criterion is that it applies 
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locally and does not rely on values at the border. In the present paper we adopt 
the picture of Aguiar et al. and our previous paper |l4j that any stable fixed 
point leads to a physical phase. Within this picture two questions arise. The 
first question is how different phases can distinguished by other properties than 
as their behaviour at the border. The second question is related to the end 
points of the disordered phase. These are second order phase transitions with 
critical indices belonging to the percolation class (l^, [l4| . To explain this some- 
what surprising behaviour we introduce besides a normal ordering parameter 
with respect to one value of the spins a second ordering on an additional group. 
By extrapolating the size of the latter to zero we will see that this ordering 
describes a bond percolation problem with a probability depending on the first 
order parameter. At the end points of the disordered phase also the percolation 
becomes critical, which explains the occurence of percolation indices. If there 
are different stable fixed points possible, they will lead to different bond prob- 
abilities. Therefore the corresponding phases can be distinguished apart from 
the boundary conditions by their different local percolation properties. 
The paper is organized in the following way. In section || we give the general 
formulae, especially for the the new order parameter. The mean field equation 
and possible chaotic solutions are discussed in section |^. In section ^ we present 
a q— state Potts model with an additional long range force allowing a physical 
interpretation of the coupling K . Section || is devoted to the appearence of 
percolation transitions in the q-state Potts model and section ^| contains our 
conclusions. 



2 General Formalism 

We consider the g-state Potts model on a Cayley tree with coordination number 
z. At each site x exists a spin a x letting q different values (f.e. a x = 1,2, ..q). 
The spins interact with their next neighbours along the bonds < xy > described 
by the following Hamiltonian 

-f3H = K Y, • (!) 

<x,y> 

Kq is proportional to the inverse temperature and the interaction strength. 
Negative Kq correspond to antiferromagnetic coupling. Partition function and 
Boltzmann distribution of a single spin can be calculated recursively with the 
methods described in [Q. A Cayley tree can be thought as z branches of 
length n connected at a site with spin er. Since the Hamiltonian (Q) can be 
written as a sum over the branches, we consider the partition function T n [a) 
for a branch where all spin summations inside the branch except the first are 
carried out. Any ordering appears as a nontrivial a dependence of T n (cr). For 
the description of this ordering we decompose the spin values into P + 1 groups 
with size qi (i=l,2,..,P+l). Each group i is characterized by a projector A 1 with 
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the property 



A* = I 1 a in gl '° Up * (2) 
' otherwise 



The group ordering of T n {a) can be written as 

P+i 

T n {<r) = V u in Al . (3) 



The disordered state corresponds to value of Ui n — u n independent of i. The 
conventional ordering with respect to a single value tr is obtained with P=l 
and A \ — b„ tS . This ordering is for the Ising case q = 2 also the most general 
one. Due to the projection property A^.Aj = 8ijA a the recursion formalism 
for P = 1 can be used also in the general case (0). The parameters Ui n obey 
a coupled system of recursion relations, as shown in the Appendix. The spin 
expectation values (A z ) lead to the magnetizations of a site. For those we use 
the linear combinations 

m i =^—(qA i -q i ) (4) 

m, = l-^-(A 4 ) . (5) 

The coefficients have been chosen such that rrii = to; = holds for a disordered 
state (A 1 ) = qi/q and mi — 1 ( m, = 1) for the ordered case where all (none) of 
the spin values a belong to group i. Among these 2(P + 2) magnetizations only 
P are linearly independent. Expressions of (A 4 ) in terms of the ratios Ui n /u\ n 
are given in the Appendix. We are mainly interested in the properties of only 
one ordering (P = 1) with respect to a group of size q = q\. As order parameter 
we use x n — 1 — U2 n /ui n . The general formula of the Appendix ( |3^ ) reduces to 
a recursion formula for x 

X n+1 =p{ Xn ){l- (l-X,^- 1 ) (6) 

with 

p(x) = (e K °-l)/(e K <>-l + qi + ( q ~ qi )(l-xy- 1 ) . (7) 

For P — 2 we treat only the case where the recursion system decouples. This 
happens for the extrapolation q^ — > 0. In addition to the recursion (|^) we obtain 
a quadratic recursion for y n = 1 — uznju\ n 

y n+1 =p{x n ){\-(l-y n ) z ~ 1 ) . (8) 

For fixed points x n = x the logistic equation (||) describes a bond percolation 
with propability p(x) fl4| . Even the order parameter y does no longer appear 
in T (since A^. — * for q% — > 0), the extrapolated value of 77T3 remains finite 
and can be related to the percolation ordering parameter mp defined by 

m P = 1 - (1 - yf (9) 
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through the relation 

m 3 =mp-(q/qi-l)mi(l-m P ) (10) 
with the magnetization belonging to the order parameter x given by 

qi {i-{i~xY) 



mi 



qi + (q-qi)(l-xY 



(11) 



(||) will be discussed in detail in the next section. The relation of percolation to 
the phase transitions of the q-state Potts model given by (|J) will be treated in 
section ||. 
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Figure 1: Iterates of the recursion ( |l4|) as function of K at go = 1-08. The 
curves give the fixed points from equ. (|17|). 



3 Pottsmodel at non integer values of q 

One of the interesting properties of (||) is the observation that it depends only 
on two combinations of the three parameters K$,q and q\. If we redefine the 
inverse temperature K by 

e K - 1 = — <e K ° - 1) (12) 

qi 
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and introduce an effective number of spin states 



qa = q/qi , (13) 

(^|) can be written as a mapping 

Xn+l = f(x n ) (14) 



with 

{e K - - {I - xf' 1 ) 
e K + (q - 1)(1 - x) z 

The recursion (^) is identical to a problem with a standard ordering parameter 
qi = 1 at temperature K with the difference that the number of states qo can 
approximate any real number qo > 1 for sufficiently large qi < q. Therefore the 
extrapolation of the q-states Potts model to noninteger values of qo is equiva- 
lent to group ordering. In the range e K — 1 > qo/z the only stable orbits of 
the recursion ( |l4| ) are fixed points. These correspond to various phases to be 
discussed in section 0. In the remainder of this section we discuss (Q) for real 
negative K and odd z. For simplicity we take z = 3. 

For qo = 1 the mapping (|lj) corresponds to the logistic equation, which exhibits 
deterministic chaos for the control parameter e K for the range 1/2 < e K < 2/3. 
For small deviations of qo from 1 we expect the mapping ( |l4| ) to have a similar 
behaviour. In figure (|l|) we show the iterates of (|l4|) for large n as function of K 
at qo = 1.08. This bifurcation diagram is similar to that of the logistic equation 
[ p"9[ . In the chaotic regions any x n is possible within the bounds 

f(l-e- K )<(l-x)<e~ K . (16) 

However, there exists a characteristic difference to the logistic equation. For the 
latter all iterates tend to oo for e~ K > 2, whereas for ( |l4| ) the chaotic regime 
ends. This can be explained by the existence of a pair of fixed points of (fL4|): 

x ± = l+ ^^ (1 - e K ± ^(1-6^-4(90-1)) (17) 
possible in the range 



K 



1| > 2^90-1 • (18) 

x+ is stable (solid line in fig ([[])) and X- (broken line) is unstable for K < 0. 
The bassin of attraction of x + is given by initial values xq > x_ . Therefore the 
chaotic regime ends where the lower bound in ( |T6| ) is equal to x- from (|l7|). 
Solving this equation for K we find the endpoint Ke of the chaotic regime 

e KE = ^(l + ^9~8q~o) ■ (19) 

If we increase q beyond 10/9 the endpoint of the chaotic regime is determined 
by equality in (iq). At the same time the pattern of the chaotic behaviour 
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Figure 2: Iterates of the recursion ([14]) as function of K at qq = 1.20. Chaotic 
regions are intercepted by stable orbits with increasing periods. The accumula- 
tion point of these orbits is given by equ (|l8|) . 



changes qualitatively. An example for qo — 1.2 is shown in figure (||). The 
regions of chaos are intercepted by stable orbits with increasing periods of r = 
2, 3, ... If we enlarge the region in between the stable orbits they are similar to 
the pattern seen in figure ([l]). The stable orbits exhibit an accumulation point 
at the end of the chaotic region given by equality in (|l8|). The period r of the 
orbit as function of K may be estimated (for sufficiently large r) by replacing 
the recursion formula by a differential equation: 



dx 
dn 



/(*) 



(20) 



The solution x{n) of (J2fj) can be found by an elementary integration. They are 
periodic with a period 



2n(l-e K )(2-q ) 



(l-e K + g ) + ^4(90 - 1) - (e K - l) 2 



(21) 



In figure (||) we compare this t{K) with the integer length of the orbits in figure 
([|). As corresponding K value we take the first bifurcation in each intervall. 
As one sees from figure (j^) the agreement between the observed periods and the 
estimate (Ell) is excellent, already at small r. For values of qo larger than 5/4 the 
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Figure 3: The continuum period (|2lj) as function of K for various values of qo. 
The data points are estimates from figures (§) , (|]) and (||) using as K value the 
locus of the first period doubling. The dashed line gives the critical curve for 
qo = 5/4. 

fixed points ( |l7| ) no longer exist. The periods as seen in figures (Q) (go = 1-255) 
and (j^) (qo = 1.275) remain finite. This is confirmed by the agreement with ( pi"] ) 
in figure @. The last chaotic region extends to K — > — oo and the number of 
regions of stable orbits decrease with increasing q . For q > q c with q c « 1.53 
there will be no transition into chaos. For 1.56 < q$ < 3 only a r = 2 orbit 
exists. For qg > 3 and negative coupling K we have only the disordered fixed 
point x = 0. 

We see that both the recursion (|l4| ) for general go and the logistic equation for 
qo = 1 have a control parameter K which decides on the details of the chaotic 
behaviour. In addition the parameter qo = q/qi decides which type of chaotic 
behaviour occurs. The various types we found are summarized in table (II]). 
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Figure 4: Iterates of the recursion (jlj) as function of K at qo — 1.255. Different 
chaotic regimes are intercepted by a finite number of stable orbits. 

4 A Generalized Potts model 

As shown in the previous section a Potts model extrapolated to non integer 
values qo is equivalent to a q-state Potts model at a rescaled temperature using 
an order parameter describing simultaneous ordering of a group of qi spins. 
The restriction to physical values of Kq imposes by ( |l2| ) a lower limit of K 1 
which excludes the chaotic behaviour discussed in section |^. Therefore still an 
extrapolation is required. To find a model, where K has a physical meaning, 
we introduce the number E of occupied links in each spin configuration: 

E = E 5 — v • ( 22 ) 

<x,y> 

A new model may be defined by allowing only spin configurations where E is 
even or odd. This corresponds to a nonlocal spin interaction. The partitioned 
sum for even E then will involve branch contributions with both even and odd 
E. It is convenient to define partitioned sums T^(a) or T°[a) for branches 
allowing only even or odd values for E in the branch. As in the previous case 
we express T^'°(a) by contributions of the subbranches T^° x {a). This system 
decouples, if we consider the sum and difference of even and odd branches given 

by 

T^{a)=T°{a) ± T» . (23) 



Figure 5: Iterates of the recursion as function of K at qo = 1.275. Different 
chaotic regimes are intercepted by a finite number of stable orbits. 

If we parametrize in the form (|32|) 



P+i 



4 ±) w = E u £ ) a; , (24) 

i 

corresponds to the unrestricted case given by ([32|) and ( |34| ) in the Ap- 
pendix. Using only one order parameter (P = 1) and introducing the ratio 



^ = %-l , (25) 



(-) 

2. ii 



we obtain the recursion formula 



-(-) - 



(i+ e ^)((i-.x(-)r i -i) 



(26) 



11+1 qi-l-e K <> + (q-q 1 )(l-x i n ) )*- 1 

Redefining the temperature by 

e K ° + l = qi (l-e K ) , (27) 

we obtain with qo = q/qx the same recursion formula as ([wj) without restriction 
on K. The chaotic region of (n4h can be reached for large qx and fixed value 
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for Kq. In this range the according recursion to x n leads always to the trivial 
stable fixed point x^ = 0. 

5 Phases of the Potts model 

The stable fixed points of recursion (^|) give in our interpretation the possible 
phases with respect to the group ordering parameter x. We do not perform 
the replacements ( |l2] ) and (|l^) , but describe the model in terms of the original 
parameters Ko,q and q\. There exist three critical points K C ,K' C and K'J. 
The disordered fixed point x = is stable at high temperatures in the range 
K' c < K < K c with K' c = ln(l - q/z) and K c = ln(l + q/(z - 2)). As the 
temperature is lowered two new fixed points given by ( |l7| ) appear for K > K". 
Only x+ with positive magnetization m\ corresponds to a stable fixed point. 
For z — 3 the value of x+ is given by equ. (jl?])- In contrast to Kq and K the 
value of K" depends on q\. For z = 3 one finds K" = ln(l + 2^q\{q — qi)). 
Therefore there exists a series of critical points depending on q\ . For Kq > K c 
the disordered phase changes into the fixed point X- ( for z = 3 given by equ. 
( |l7|) ) with negative magnetization mi, and for Kq < K' c we encounter an 
antiferromagnetic coupling. This transition corresponds to the first bifurcation 
discussed in section ||. If more than one stable fixed point exists, the value at 
the boundary decide about the adopted phase. The transitions at K" is of first 
order. For the following it does not matter wether this transition apears at K" 
or a higher value as advocated in |J from equality of the free energies. It is 
only important that the transitions at K c and K' c are genuine transitions and do 
not correspond to spinodal points. For the critical indices at K c we find values 
belonging to the percolation class ((3 = 7 = 1). The same is true for the average 
magnetization also at K' c . The staggered magnetization equal to the difference 
of magnetizations at adjacent sites exhibits at K' c indices of the mean field class 



q range 


end point exp(X) of 
chaotic behaviour 


repetition of 
deterministic chaos 


maximum 
orbit length 


1 < qo < 10/9 


1/4(1 + V9 - 8<fo) 


1 


00 


10/9 < q < 5/4 


1 - 2^90 - 1 


00 


00 


5/4 <q <q c 





2^(2-go)/(go + l)V4<Zo -5 


00 


qa > q c 







finite 


1.56 < q < 3 
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Table 1: Characteristics of chaotic behaviour in different intervalls of go- The 
value of q c above which no orbit of infinite length occurs is given numerically 
by q c = 1.53. 
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(2/3 = 7 = 1)- The Ising case qi — q/2 is exceptional. Since the amplitude of 
the average magnetization vanishes due to the Ising symmetry on an A£?-lattice 
(m(x+) +m(x-) — 0) the percolation transitions dissappear, and one observes 
mean field indices at both critical points K' c and K" = K c . Note that a Ising 
type behaviour does not only occur for q = 2qi = 2 but also for any even q with 
the choice qi = q/2. This reflects the fact that the recursion (||) can be written 
in the form ( |l4|) which depends only on the ratio q = qjq\. 
The appearance of percolation indices for q x ^ q/2 is somewhat surprising. To 
understand the mechanism we consider the recursion (jfy for the order parameter 
x together with an extrapolated second order parameter y n given by (||). The 
recursion for the latter reads at a fixed point for x as 

y n+ i =p(x){l-{l-y n y- 1 ) (28) 

with p(x) given by (0). As shown in [Q the recursion (|2^) is equivalent to bond 
percolation with a probability \p{x) | . A non zero fixed point value y leads to a 
finite probability mp for a site belonging to the infinite cluster 

mp = 1 - (1 - yf . (29) 

The mapping ( p8| ) has the fixed points y — and y — x. The stable fixed point 
must satisfy 

\p(x)(z-l)(l- y y- 2 \<l . (30) 

Both the disordered fixed point x = in the range K' c < Kq < K c and X- with 
negative magnetization for Kq > K c lead to the nonpercolating state y = 0. In 
contrast the fixed point x+ requires y — x + . Therefore the two phases possible 
for Kq > K" can be distinguished locally by their percolating properties. The 
phase with mi > exhibits an infinite percolation cluster, whereas the phase 
with mi < is in the disordered percolation state. Using the criterion of ex- 
istence of an infinite cluster the phase on the Cayley tree is uniquely defined 
and can be carried over to the infinite Bethe lattice, since one does no longer 
need the values at the boundary. For Kq < K' c the recursion (^) exhibits an 
r = 2 orbit corresponding to antiferromagnetic ordering. By a straightforward 
calculation one can show that the recursion (||) has a stable fixed point y = 0, 
although p alternates between p(x±) in subsequent iterations. 
If Kq approaches the end points K c or K' c of the disordered phase, \p(x)\ ap- 
proaches the critical value l/(z — 1). At these points all moments of the cluster 
size distribution diverge. These moments can be represented by derivatives of 
the magnetization mp with respect to an external field. Therefore magneti- 
zation and susceptibility will be singular with critical indices belonging to the 
percolation class. Since the various magnetizations are not independent due to 
relation ([l0]) these divergences will show up also in other magnetizations than 
mp. 
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6 Conclusions 



We introduced a new type of ordering in the q— state Potts model where a spin 
belongs to one or more groups of values instead of the conventional ordering 
with respect to one value. The recursion formalism for the latter can be also 
applied for the new ordering. For only one group of size q\ the recursions can 
be cast into the same form as for gi = 1 with the exception that the effective 
number of states is given by qo — q/qi , which can approximate any real number 
larger than 1. This allows an extrapolation of the q— state Potts model to non- 
integer values of q without losing the physical interpretation. For qo close to 1 
the recursion exhibits deterministic chaos similar to the logistic equation. The 
control parameter is essentially the temperature K^ 1 . Increasing go we observe 
as function of K repetitions of deterministic chaos separated by stable orbits of 
length t = 2, 3, 4, • • •. Number of the orbits and their length can be estimated 
by a continous length approximation. Most of this chaotic behaviour occurs 
at K(Ko) values, which cannot be reached with real values of Kq. To avoid 
this extrapolation to unphysical values of the temperature one can introduce 
a generalized Potts model, where possible spin configurations are restricted by 
the condition that the number of links connecting sites with equal spin values 
must be even or odd. In |hJ the Potts model has been extrapolated to real 
values of q using the analytic expression for the free energy. For q < 2 negative 
values for the specific heat have been found. This unphysical behaviour can be 
avoided if one uses the new order parameter referring to a group of q\ spins and 
interprets the ratio go = q/qi as effective number of spins. 
The phase structure of the g— state Potts model with ordering with respect to 
one group of size gi are qualitatively similar to the conventional gi = 1 case. 
Only for even g the Ising case gi = g/2 behaves differently For gi ^ g/2 the 
transitions at the endpoints of the disordered phase are of second order with 
percolation class critical indices. To understand their occurence we added a 
second group of ordering with size extrapolated to zero. This does not change 
the recursions for the first group. In addition one obtains a logistic equation 
for the second order parameter, which leads to indices of the percolation class. 
Since the control parameter in the logistic equation depends on the phase, dif- 
ferent phases can be distinguished by their percolation property. The phases 
with positive magnetization exhibit an infinite cluster with finite probability, 
whereas negative or zero magnetizations < have zero probability. This crite- 
rion allows a local distinction between two possible phases. Since one does not 
have to resort to boundary conditions, this criterion can be applied also on the 
infinite Bethe lattice. 
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stable orbits 10 
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A Recursion on a Cayley tree 

The partitioned sum for a Cayley tree can decomposed into z branches T n (a) 
of length n connected at a central site with spin a. Each branch can be related 
according to figure (|^) to z — 1 branches of length n — 1 by 

T n (?) = X>xp(Xo<W) CTn-iOr')) 2 " 1 (31) 
Parametrizing T n (a) with 

T n (a) = J2 u i,nK (32) 

i 

the products in ( |3l|) can be evaluated using the relation 

A^Ai = (33) 
Comparing coefficients on both sides we find 

Ui, n = (e K ° - 1) (li^.i) 2 - 1 + J2<U Kn-i)*" 1 ( 34 ) 

3 

Spin expectation values can depend only on the ratios of Ui <n . For these we 
obtain the recursion 



(e Ko ~ 1) 1 - l) 



"1," e K ° ~ 1 + J2j 1j («j,n-l/«l,n-l) 



2-1 



(35) 
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Figure 6: Recursion for z = 4 



Setting U2.n/u\ t n — 1 — x n and U3 lV ,/v,i in = 1 — y n we obtain at (73 = the 
recursions (|) and (|) in the case P=2. The Boltzmann weight for the central 
spin at a fixed point = Uj/iti is given by 

w(cr) = 7V(T(ct)) z (36) 

where the normalization constant follows from iu(<t) = 1. With the distribu- 
tion @ we get for (A 1 ) 

(A*) = ifr )Z . . (37) 

<& + £ 3 >i<&fa) a 

If we use in the case P = 2 as independent magnetizations mi and 77X3, we find 

mi = «i(l-(l-a:)*)/iV (38) 

and 

m 3 = 1 - q{l - v) z /N (39) 
with the normalization factor N = q± + (q — qi)(l — x) z . 
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